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We study erratically moving spatial structures that are found in a driven interface in a random 
medium at the depinning threshold. We introduce a bond-disordered variant of the Sneppen model 
and study the effect of extremal dynamics on the morphology of the interface. We find evidence 
for the formation of a structure which moves along with the growth site. The time average of the 
structure, which is defined with respect to the active spot of growth, defines an activity-centered 
pattern. Extensive Monte Carlo simulations show that the pattern has a tail which decays slowly, as 
a power law. To understand this sort of pattern formation, we write down an approximate integral 
equation involving the local interface dynamics and long-ranged jumps of the growth spot. We 
clarify the nature of the approximation by considering a model for which the integral equation is 
exactly derivable from an extended master equation. Improvements to the equation are considered by 
adding a second coupled equation which provides a self-consistent description. The pattern, which 
defines a one-point correlation function, is shown to have a strong effect on ordinary space-fixed 
two-point correlation functions. Finally we present evidence that this sort of pattern formation is 
not confined to the interface problem, but is generic to situations in which the activity at succesive 
time steps is correlated, as for instance in several other extremal models. We present numerical 
results for activity-centered patterns in the Bak-Sneppen model of evolution and the Zaitsev model 
of low-temperature creep. 

PACS numbers; 05.40. -fj, 47.54.-(-r, 68.10.Gw, 47.55.Mh 



I. INTRODUCTION 

Driven interfaces in random media present several fea- 
tm'es of interest, both with regard to the morphology 
of the moving interface as well as the dynamics of the 
growth process. Experiments have been performed on 
several sorts of systems, ranging from fluid flow in porous 
media [0-|5j to propagation of burning fronts These 
indicate that the disorder in the medium affects the prop- 
erties of the interface in a crucial way. In particular, the 
large-distance scaling properties differ considerably from 
those of interfaces in uniform media. In both theoret- 
ical and experimental investigations, it is customary to 
characterize the spatial structure of the interface by its 
roughness. The main point of this paper is to show that 
there is sometimes an unusual sort of pattern formation 
in the system, which results in the interface acquiring 
a time-averaged shape. In such situations, this pattern 
provides an alternative characterization of interface mor- 
phology. 

A customary measure of the roughness is provided by 
the exponent a, defined by ly ^ L", where W is the 
root mean squared width of the interface and L is the 
size of the system. The experiments mentioned above, 
and others similar to these, report an anomalously large 
value of a — large compared to the predictions of ex- 
isting theories for interface growth in nonrandom media 
1^ . It is recognised that the quenched nature of the dis- 
order in the medium is responsible for this difference in 



scaling properties of the interface. Unlike thermal noise, 
which varies rapidly, a portion of the interface subject 
to a quenched-noise environment continues to experience 
the same forces until the growth process takes it forward 
to a new region. The pinning effect of quenched noise 
has a strong effect on the large scale properties of the 
interface. 

Several theoretical models have been put forward to 
try to account for the effect of quenched disorder on 
the properties of the interface; j|] gives an account of 
some of the early work and the relationship to other 
problems involving pinning, while [p^ is a a recent re- 
view. Amongst various proposals put forward to ex- 
plain anomalous roughening are theories based on con- 
tinuum equations with quenched disorder pd|-[lq|, the 
inclusion of noise with power-law amplitude |17| , |18| or 
long-ranged correlation ||l9H21|, power-law distributions 
of pinning-center strengths ]22|| , as well as a class of 
models with microscopic rules based on directed perco- 
lation |^,^,|4| and models which relate the large-scale 
structure to the wetting properties of the invading fluid 
p6t . A number of these models base the explanation 
of anomalous roughening on the phenomenon of critical 
depinning, which is relevant to an interface just at the 
threshold of motion. In the opposite limit of large ve- 
locity, the interface encounters the disorder at any site 
only for a short time, suggesting that the quenched na- 
ture of disorder is not important in this limit, and the 
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FIG. 1. A schematic picture of the activity-centered pat- 
tern in an untilted {m — 0) interface showing (a) the height 
gradient pattern as defined in Eq. 1, (b) the interface profile 
obtained by integrating the pattern in (a). The discontinu- 
ity at the origin in (a) leads to a cusp at the origin in (b), 
implying that the active site is most likely at the peak. 

interface behaves much as in a nonrandom medium . 
However this is not true at low velocities of the inter- 
face near the depinning threshold. In particular, the 
limit of zero velocity is thought to be a dynamical critical 
point 1^ where the scaling properties of the interface are 
strongly affected by the disorder. 



In a certain class of models, the quenched disorder en- 
ters as barriers of random strengths which impede inter- 
face motion. The formation of infinitely long directed 
percolating paths of these barriers is of special signifi- 
cance, as such paths can block the entire interface effec- 
tively The model proposed by Sneppen HI] 
involves "extremal" dynamics: At each time step, the 
interface advances only along the weakest of the barri- 
ers. Extremal models were first introduced much earlier, 
in the context of invasion percolation for two-phase fluid 
flow in porous media, with a non- wetting fluid displacing 
a wetting one [27|-|2g|]. The predictions of this invasion 
percolation model were borne out by experiments 
Dynamical correlation functions involving the center of 
activity in invasion percolation obey scaling ]3l]j3^ , and 
the process defines a self-organised critical phenomenon 
p2t ; the interface organises itself to align along critical 
paths bordered by large barriers, without the necessity 
of tuning any external parameters. The model proposed 
by Sneppen p^ ] is a modification of the invasion perco- 
lation model, incorporating surface tension effects which 
prevent very strong local convolutions of the interface, re- 
sulting in a self-affine, rather than self-similar, geometry. 
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FIG. 2. A schematic picture of the activity-centered pat- 
tern in a tilted interface (m 7^ 0), in (a) height gradients (b) 
the interface profile obtained by integrating the pattern in 
(a). As can be seen in (b), there is a larger than nominal 
slope near the active site. 

An extremal model close to the Sneppen model was 
shown to result js^] from a model of wetting-fluid in- 
vasion, by considering an interface advancing by merging 
meniscus arcs between adjacent pairs of pinning centers. 

An interesting feature of extremal dynamics is that it 
induces strong spatial correlation in sites at which growth 
occurs p^-|36t at successive time steps. In this paper 
we introduce and study a variant of the Sneppen model 
and show that the interface develops an interesting time- 
averaged structure as a result of correlations. The defin- 
ing equation for the structure is 

'f{r) = ^[(Wh(r + R{t)))-m] (1) 

where R{t) is the position of the active site at time t, 
h{r') denotes the height at the site r', (...) is a time aver- 
age in the steady state, and m is the overall slope of the 
interface. The unusual point is that this structure ^'(r) 
is not fixed in space, but moves with its center always at 
the growth site, which itself follows an erratic path. The 
moving origin is crucial to the definition, as time averages 
performed at a fixed point in space reveal no structure 
at all. This average structure defines a pattern, and we 
study its formation numerically and analytically. 

Figure la depicts schematically the activity-centered 
pattern in height gradients with respect to the moving 
origin, for the untilted (to — 0) interface. As can be 
seen, the interface develops an overall shape which is de- 
scribed by the height-gradient profile ^(r) . The tail of 
the pattern falls slowly at large distances, as a power 
law. The corresponding height profile of the interface 
h{r) with respect to the active site is represented in Fig. 
lb. The nature of the pattern is sensitive to tilt, and 
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Figure 2 shows the height-gradient pattern and height 
profile for a tilted interface. 

We propose that this pattern is a simple way of char- 
acterising a new aspect of the morphology of the inter- 
face. Traditional ways of characterising the morphology 
involve, as has already been mentioned, determining the 
roughness exponent a. However such a definition does 
not hone in on the overall shape of the interface. In sit- 
uations such as the one considered in this paper when 
the interface does develop a nontrivial structure, the pat- 
tern is a useful quantitative characterization. An impor- 
tant point about the time-averaged pattern is that Eq. 
1 defines a one-point correlation function. As such, it 
would be expected to strongly influence the properties 
of customary two-point correlation functions. We ver- 
ify this by numerically studying two-point correlations. 
Further we find that this sort of pattern formation is not 
restricted only to the Sneppen model, but also occurs 
in other extremal models, such as the Zaitsev model for 
low-temperature creep and the Bak-Sneppen model 
p8| of biological evolution, albeit in other quantities. 

The plan of this paper is as follows. In Section II we 
introduce our model and discuss the connection with the 
problem of directed percolation, well established from 
earlier studies. In Section III, we discuss the correla- 
tions in the location of successive growth sites, a concept 
central to this paper because of its connection with pat- 
tern formation. In Sections IV and V we define the pat- 
tern and present numerical results as well as an integral 
equation which provides an understanding of this sort 
of pattern formation. We define a model for which the 
equation is exact and discuss how the approximation can 
be improved. In Section VI, we discuss the issue of tem- 
poral correlations with a view to seeing how they affect 
the pattern. In Section VII, we present our results for an 
ordinary two-point correlation function in our model and 
show that activity-centered pattern formation has to be 
taken into account in order to understand some features 
in it. Section VIII deals with pattern formation in other 
extremal models and we conclude with a summary of our 
results in Section IX. 



II. EXTREMAL MODEL OF INTERFACE 
DEPINNING 

The extremal-model description of fluid-fluid interfaces 
in porous media is valid when the wetting is dominated 
by capillary forces, and thermal fluctuations are not im- 
portant. In the extremal model, the interface advances 
along the weakest barrier just ahead of it. The ap- 
pealing feature of the model is that it is self-organised 
critical; the dynamics, which involves searches for the 
global minimum at every step, automatically tunes the 
interface to a critical state at the depinning transition, 
without the necessity of fixing any external parameter. 




FIG. 3. The Extremal Bond Model. The tilted interface 
//' advances along the extremal perimeter bond A and locally 
readjusts to align along the dashed line. At the next instant, 
the activity moves from A to A'. The corresponding config- 
uration and local moves for the particle-hole model are also 
shown. SS' is a stopper, whose perimeter is fully occupied by 
diodes. 

In the model proposed by Sneppen ||2^ , the random 
medium is modelled by a square lattice in which the sites 
are assigned random numbers / G [0,1]. The random 
numbers could signify, for instance, the pore sizes in a 
porous medium. The interface is a directed path on this 
lattice and grows only at that perimeter bond with the 
smallest value of the random number; after every such 
move, a local rearrangement process ]39| ensures the ab- 
sence of very large slopes. 

Extremal dynamics has also been proposed to describe 
very different situations — for instance, the phenomenon 



37y33f| , crack prop- 
38^. We will see in 



of low-temperature dislocation creep 
agation and biological evolution 
Section VIII that the sort of pattern formation we find 
in the interface model occurs in these models as well. 



A. The Extremal Bond Model 

We study a modified version of the Sneppen model in 
this paper. In this version, hereafter referred to as the 
Extremal Bond model (EBM), the interface is taken to 
be a directed path on a square lattice (Fig. 3), with 
tilted cylindrical boundary conditions pd| which ensure 
that the mean slope is preserved. To every bond k on the 
lattice is pre-assigned a fixed random number fk drawn 
from the interval [0,1]. The interface grows at that bond, 
in front of the interface, which carries the smallest ran- 
dom number. The local growth rules are the following; if 
the chosen minimal bond has a positive (negative) slope, 
the sequence of links with negative (positive) slope just 
below (on the left) also advances, as illustrated in Fig. 
3. This preserves the length of the interface as would 
happen in situations with very high surface tension. The 
local dynamics of interface adjustment is similar to that 
for the low-noise Toom interface in |^^; the models are 
different in that the growth site is picked by the extremal 
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rule in our case, while it is picked stochastically in the 
Toom interface model. 

The EBM differs from the Sneppen model in that the 
length of the interface is a strict constant of the motion 
and the fkS are associated with bonds rather than sites. 
While this modification does not change the values of any 
of the large- distance scaling properties of the Sneppen 
model, it has a few advantages. The interface aligns along 
directed spanning paths in a percolation problem, just 
as in the Sneppen model. In our case, the correspond- 
ing percolation problem is the diode-resistor percolation 
problem [Q. On the square lattice, it is dual to the di- 
rected bond percolation problem |Q which is relatively 
well-studied. Another advantage is that the problem of 
interface growth in the EBM is conceptually simplified 
by the existence of a known one-to-one correspondence 
between the growing interface and a system of hard-core 
particles moving on a ring. The two-dimensional prob- 
lem hence reduces to an effectively one-dimensional one. 
This also facilitates the numerics. The correspondence 
between the interface and the hard-core particles is de- 
tailed below. 

Positive slope links of the interface are represented 
by particles {rij — 1) and negative-slope links by holes 
{fij — 0); see Fig. 3. The difference in height of the 
interface between sites ji and j2 is given by hj^ — hj^ = 
X^jLji — !)■ In front of each link of the interface is 
a bond with a random number / assigned to it. Corre- 
spondingly, the site j with the particle (hole) represent- 
ing this link carries a random number fj . Just as for the 
interface, at each time step, activity is initiated at the 
site with the minimum fj. The update rules for the in- 
terface translate to the following dynamics for particles 
and holes. If the site with minimum fj contains a par- 
ticle (hole), it exchanges with the first hole (particle) to 
the left (right). All sites hopped over, including the two 
which exchange the particle and hole, are refreshed by 
assigning a new set of fj's. This corresponds to the fact 
that the updated portion of the interface moves ahead 
and meets a fresh set of /'s on the square lattice. Be- 
cause the number of positive slope links (and hence also 
the number of negative slope links) is conserved for the 
interface, in the particle-hole terminology, this implies 
that the number of particles is conserved. Hence we can 
define a density p for particles on the one-dimensional lat- 
tice. This density determines the mean slope m — 2p ~ 1 
of the interface. An untilted interface corresponds to 
half-filling. The reference direction for determining tilt 
is the easy direction of directed bond percolation on the 
square lattice, which is along the 45° line. Tilt refers 
to any density away from 0.5, which implies a slope dif- 
ferent from 45°. The interface advances in a direction 
perpendicular to the direction of tilt and this translates 
to a nonzero current of particles on the ring. 



B. Connection to Diode-Resistor Percolation and 
Directed Percolation 



Extremal models of interface depinning make use of a 
correpondence to the Diode Resistor Percolation (DRP) 
and Directed Percolation (DP) problems to predict vari- 
ous properties of the interface. In view of this, it is useful 
to recall some facts about the DP and DRP processes. 

In the directed percolation problem, bonds on a lattice 
are occupied with probability p. At some critical value pc 
an infinite directed path of occupied bonds (in which ev- 
ery step is taken rightward or upward) first forms along a 
definite direction; on a 2 — c? square lattice, this is along 
the 45° direction. For p > pc, The network of these 
infinite paths forms an infinite connected cluster. For 
directed bond percolation on a square lattice the value 
of Pc is known to be ~ 0.6446 |4^. There are two dis- 
tinct correlation lengths, along the easy direction and 

transverse to it, both of which diverge as p — > 
Cj| {p — Pc)~''^^ ,£,± ~ {p — Pc)~'^^ respectively. The 
values of these exponents are known to be ~ 1.733 
and lyjL ^ 1-097 @. 

Suppose we have a single source point, and we ask 
which portion of the plane can be reached from it via 
occupied directed bonds. For p > pc this region is con- 
tained within a cone with opening angle 2<j) = arctan(m) 
where m is the slope of the edge of the cone with respect 
to the 45° direction; the opening angle depends on p. 
This relation can be inverted to find the critical proba- 
bility Pcim) viz. the probability at which a connection 
first appears along the direction with slope m ^ 0. Cor- 
relation lengths along and perpendicular to this direction 
have exponents — 1 and iy±_ — 0.5 |Q. We refer to the 
direction along 45° as untilted (to = 0); any other slope 
is referred to as tilted. 

In the diode-resistor percolation problem, every bond 
is occupied by a "diode" (a one-way connection) with a 
probability p or a "resistor" (a two-way connection) with 
a probability 1 — p. On a square lattice, the diodes all 
point up or right. Let us ask which regions of the plane 
are connected to a given source point. If p = 0, a source 
point can reach the entire quadrant of which it is the left 
corner. As p decreases from 1 to Pc, the opening angle 
4>' of the connected region increases from 7r/2 to tt; be- 
yond this, the entire plane can be reached from the source 
point. The edge of the connected region is bordered by 
diodes pointing rightward and upward, which prevent it 
from spreading leftward and downward (Fig. 4). 

On the 2 — d square lattice, DP and DRP are 
dual to each other The dual to a DRP con- 

figuration is constructed using the following rules. 
A diode in the DRP lattice is crossed by a diode 
in the dual lattice, whereas a resistor is crossed 
by an insulator (no connection) in the dual lattice. 
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FIG. 4. The boundary of the region reachable from the 
origin in diode-resistor percolation. For a nonzero fraction of 
resistors, the opening angle of the region is larger than 90° as 
shown. The dashed line is an infinite path in the correspond- 
ing directed percolation problem on the dual lattice. 

Thus we recover the directed percolation problem on the 
dual lattice. The opening angle of the cones in the two 
problems are related hy (j) + 4>' — n . 

In the EBM, the random medium is modelled by con- 
sidering a square lattice with every bond assigned a ran- 
dom number fk drawn from the interval [0,1]. For a 
certain trial value /*, imagine occupying all bonds with 
/ < /* by resistors, and the rest with diodes. We thereby 
generate a DRP configuration with p = 1 — /*. When /* 
takes on the value I — pc, an infinite connected path of 
diodes is formed. Such a path is called a 'stopper', and 
is significant for the dynamics of the EBM, as a moving 
interface with no overall tilt will align with such stop- 
pers from time to time When the interface aligns 
along a stopper, all the bonds in front of it are larger 
than /c = 1 — Pc (Fig. 5a). Similarly a tilted interface 
with slope m aligns along the edge of the cone with the 
same slope and all the bonds in front of it are expected 
to have a value larger than fc{m) = 1 — pc{m). Since 
p,{m) > p,{0) for m ^ 0, fc{m) < /,(0) (Fig. 5b). 
Even when the interface is evolving between two stop- 
pers, only a small fraction of its overall length actually is 
in between; the rest is still aligned with a stopper. The 
non-aligned fraction is expected to vanish in the ther- 
modynamic limit. These expectations are confirmed by 
numerical studies of the EBM. As can be seen from the 
figures, the bonds in front of the interface are all mostly 
larger than a threshold value. 

Consider an interface of slope m aligned 
along a critical DRP path of the same slope. 
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FIG. 5. The steady state probability P{f) of having a 
bond with the value of the random number equal to /, in 
front of the (a) untilted (b) tilted interface with p = 0.75, in 
the Extremal Bond Model. There is a tilt dependent thresh- 
old below which the probabilty of finding a bond with that 
value is very low. A system of size 1000 was averaged over 
lO'^ configurations 

It then moves forward by puncturing the path at the site 
with the least value of /, which for an infinite system 
is exactly 1 — pc{m). On piercing through, a portion of 
the interface grows and fills out a loop of the infinite 
cluster while the rest of it remains pinned. However the 
interface motion within the loop is far from uniform. 
Just as the critical cluster at pdjn) impedes the growth 
of the interface on length scales of the order of the size 
of the system, near-critical clusters impede its motion at 
length scales of the order of but smaller than the loop 
size. One can think of these clusters as forming a finer 
network of connections within the network formed by the 
critical cluster at Pc{m). While the interface is filling out 
a loop of the critical cluster, it encounters this finer mesh 
and as a result its motion is impeded temporarily. In 
what follows, we refer to these near-critical connections 
as 'sub-stoppers'. A substopper can be characterized by 
the lowest value of fk on it, say f^^, and also by the 
typical length scale / over which it provides for effec- 
tive pinning of the interface. These are related through 
ir^-/c(m)|-'^ii ~L 
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FIG. 6. The location of the active site as a function of time 
for (a) untilted (b) tilted interface. The region marked (1) is 
a typical example of a window in time in which the activity 
is localised in a region. In (2) the activity ranges through the 
whole system, corresponding to an instance when the inter- 
face has been pinned by a critical cluster at pc{m). Though 
the directionality of the active-site motion is evident in (b), 
this does not induce a net drift as mentioned in the text. The 
data displayed above is for L = 1000. 



III. CORRELATIONS IN THE ACTIVE-SITE 
MOTION 

The above description of the evolution of the interface, 
contained as it is by networks of substoppers and stop- 
pers, makes it clear that there are strong correlations 
between successive points of growth or forward motion. 
These correlations extend from small length scales up to 
scales of the order of the system size. Figure 6 shows the 
plot of the location of the active site for 10,000 time steps 
for both the tilted and untilted cases. It can be seen that 
there are jumps on all scales in the active site position. 
The figure corroborates the description of interface mo- 
tion given above. Region (1) is a typical instance of the 
interface filling out a loop of size I. It shows that there 
are jumps of all sizes upto a length I, bearing out the 
substopper picture. Region (2) on the other hand marks 
an instance when the interface has aligned along a stop- 
per Pcim) and hence there are jumps of all sizes upto the 
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FIG. 7. Monte Carlo results for the probability distribu- 
tion of the jump of the active site for three different densities, 
p — 0.5 (plus sign), p = 0.75 (circles) and p = 0.84375 (tri- 
angles). If p 7^ 0.5, p{l) is not a symmetric function and p(l) 
and p{—l) (both of which are shown in the figure in the two 
p 7^ 0.5 cases) do not coincide as they do for the symmet- 
ric case. While p{l) > p{—l) for small I, the curves cross so 
that the situation is reversed for large jumps. The two curves 
asymptotically coincide with an asymptotic slope that differs 
from that for p — 0.5 We used L= 65536 and averaged over 
3.10^ configurations. 

system size. 

More quantitatively, a measure of this long-ranged mo- 
tion of the active site is the probability distribution p{l) 
that two consecutive locations of the active site are a 
distance / apart. Figure 7 shows p{l) for both the tilted 
and the untilted interface. In both cases, p{l) decays as 
a power for large I: p{l) ^ |^|^". 

In the untilted {p = 1/2) case, p{l) is a purely sym- 
metric function because of r ^ —r symmetry. We find 
TT — 2.25 ± 0.05, which compares well with earlier deter- 
mined values of vr for the Sneppen model p^-|36[|. 

In the tilted case {p ^ 1/2), p{l) is not a symmet- 
ric function (Fig. 7). As can be seen from the figure, 
there is a larger number of small jumps to the right, 
but more jumps of large magnitude to the left. It is 
convenient to separately analyse the even and odd parts 
p± = {p{l) ±p(— Z))/2 in order to find the exponents. We 
find that the even part, p+{l), decays asymptotically as 
p+{l) ~ |/|-''+ with7r+ = 2.00±0.02. The odd part p_(/) 
changes sign (as implied by the crossing of the curves 
in Fig. 7) and asymptotically follows p-{l) ^ 
with TT- = 2.49 ± 0.06. We verified that the values of 
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FIG. 8. (a) The function p{l'\l) as a function of I is 
logarithmically binned in powers of two and is maximum for 
the curve with the largest flat stretch. As can be seen in the 
figure, the flat stretches of the large-Z curves coincide, peeling 
off at a value which is I dependent, (b) The function p{l'\l) 
plotted as a function of I' /I. The curves coincide for large I 
indicating that the conditional probability is just a function 
of the ratio I' /I. 

7r+ and 7r_ are the same for various p ^ 1/2. An inter- 
esting aspect of the function p{l) for the tilted interface 
is that J lp{l)dl vanishes in the thermodynamic limit. In 
terms of the active-site motion, this implies that though, 
for the tilted interface, the short-ranged jumps of the 
active site are mostly along the tilt, there are enough 
long-ranged jumps in the opposite direction to balance 
this. The cancellation gets better as the thermodynamic 
limit is approached. The tilt of the interface, therefore 
does not induce a net drift in the active site motion. 

Because of the close connection of the interface grovirth 
problem to DP, several exponents associated with the 
growth and structure of the interface are thought to be 
related to the DP exponents v\\ and v±_ by scaling rela- 
tions 35,3^]. The exponent tt in the untilted interface 
as well as the exponent 7r+ in the tilted case have been 
argued to be related to the DP exponents by the scaling 
relation tt = 1 -t- (1 -I- v±_)/v\\ |^^. Besides this, for the 
tilted case, the exponents can also be obtained exactly 

As mentioned earlier, the function p{l) does not con- 



tain all the information about the active-site motion. 
This is due to the presence of temporal correlations in 
the interface depinning model: the jump of the active 
site at the present instant is strongly correlated to the 
jumps before. As a result, what is needed is the full 
distribtion function ^2!^!) which is the condi- 
tional probability distribution that a jump of length I 
occurs ai t = n given that a jump of length li occurred 
at i = 1, a jump of length I2 occured at i = 2 and so 
on. The probability distribution p[l) is obtained by in- 
tegrating out the other variables li,l2----ln-i- In order 
to assess the importance of these temporal correlations, 
we numerically determined the conditional jump proba- 
bility distribution p{l'\l), the conditional probabilty that 
a jump of length /' occurs given that at the previous in- 
stant, the active-spot jumped a distance I. Figure 8a 
presents our numerical measurement of this function for 
the EBM. After the initial short-ranged decay, which is 
a result of the local back- wetting move of the dynamics, 
the function is flat over a considerable range, and then 
decays as a power. As evident in the figure, the range of 
the flat region depends on / but the value of the probabil- 
ity in the flat region becomes independent of / for large /. 
In Figure 8b we show the scaling plot of p{l'\l) obtained 
by plotting the curves as a function of V /I. The function 
evidently approaches a scaling form for large I and I' . The 
scaling function is flat over a region and decays beyond 
as a power ~ {I' /l)^^ . We expect tt' = 7r(~ 2.25) to 
hold as I' /I 00. Although the measured value of tt' in 
the range shown is larger (~ 2.9), the bending apparent 
in the lower right portion of the curve is consistent with 
an approach to the value tt. 

Qualitatively, the behaviour of this function may be 
understood thus. If the active site jumps a distance I 
at the previous instant, one can think of the interface as 
pinned by a DP cluster with loops of average linear di- 
mension ^ I. Most of the interface would then be pinned 
while a portion of it fills out a loop of linear dimension 
^ I. This would imply that on an average, any jump 
smaller than / is equally likely. On length scales larger 
than I the motion is like the original problem and hence 
the jump probability decays as a power. This line of 
argument would imply that p{l'\l) should be a scaling 
function of I' /I. Figure 8b bears out this expectation. 

Another equivalent way of understanding the function 
p{l'\l) is using the 'backward-avalanche' technique intro- 
duced in [Q. A backward-avalanche is defined as fol- 
lows. If at time s + S a random number fk is picked as the 
minimum, the magnitude of the backward avalanche is S 
if at time s the random number picked was larger than 
fk but for time s + r for r < S the value of the random 
number picked was smaller. That is, one goes back in 
time to the first instant when the random number picked 
is larger than the present one and this interval of time 
is the magnitude of the backward-avalanche initiated at 
the present instant. If the active site hopped a distance 
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FIG. 9. The distribution F{s\l) of durations of avalanches 
for jumps as defined in the text. For small I, F{s\l) decays 
exponentially (filled triangles) as expected, since it is very un- 
likely that a larger jump does not immediately occur. As I 
increases, the cutoff for the decay increases (empty and filled 
circles) till for large enough I, the curves coincide, decaying 
as a power ~ s~^. 

I in the previous instant, this would imply that the size 
of the spatial region affected by the backward avalanche 
initiated at the previous instant is /. At the next time 
step any of three possibilities can occur. If the random 
number picked is smaller than the previous one, then the 
active site necessarily lies in the local region which has 
just been affected by the move. If the random number 
picked is larger, then either the random number picked 
belongs to the same backward-avalanche as the previ- 
ous one, or it belongs to a bigger one (which therefore 
encompasses the previous one). In the former case, the 
active site is equally likely to hop anywhere in the region 
affected by the present backward-avalanche ( [^,Q) 
and hence the function is flat upto a length ~ L In the 
latter case, active-site hops are in general larger, since 
the encompassing backward avalanche is larger. The 
probability of having the active site hop a distance li > I 
is then related to the probability of finding a backward 
avalanche of this spatial extent, and hence decays as a 
power. 

Another indication of the correlation in the jumps of 
the active site is the following quantity which is the ana- 
log of an avalanche, for jumps (Fig. 9). Let the avalanche 
be initiated at an instant s by a jump of magnitude / of 
the active spot. The "jump-avalanche" lasts as long as 
consecutive jumps in the active site are less than I. That 
is, the avalanche is of duration S if at time s the active 
site hopped a length I, and for S consecutive instants af- 
ter that, the hops in the active site are smaller than I 
till at the (5 + l)th instant, it hopped a distance greater 
than I. The analog of this quantity for random numbers 
was first defined for models in the non-wetting invasion 



percolation regime [ pT|,p2[ and later also measured for 
the Sneppen model P5| , p5| , ^ . The jump-avalanche is a 
strong indication of the correlation in successive jumps of 
the active site. As Fig. 9 indicates, the function falls off 
as a slow power. This is an indication of the long-term 
memory in the system, as the function would decay ex- 
ponentially if successive jumps had not been correlated. 



IV. PATTERN FORMATION IN THE 
INTERFACE DEPINNING MODEL 

A. Definition 

As discussed in the previous section, the motion of the 
active spot in this model is very correlated. One can then 
ask whether this motion reorganizes the shape of the in- 
terface (or equivalently the arrangements of particles and 
holes in the particle model corresponding to the EBM) in 
any specific way. Below we will show that the correlated 
motion of the growth-spot in this model leads to the for- 
mation of a pattern in height gradients (or equivalently 
the density of particles) of the interface. 

The defining Eq. Q for height gradients can be written 
in terms of densities as 



*(r) = {n{r + R{t))) - p. 



(2) 



Here n(r) is the density at site r and the angular brackets 
denote a time average in the steady state. 

This pattern is linked to the formation of a structure 
which is very different from normal space-fixed patterns 
because it is referred to an origin R{t) which is moving 
around. In fact it can be discerned only when the origin 
moves. An ordinary space fixed average is translationally 
invariant and satisfies {n{r)) — p for any site r. Since the 
pattern is centered around the site of activity, we will re- 
fer to it as the "Activity-centered pattern" (ACP). 



B. Numerical results 

We studied ^'(r) by Monte Carlo simulation in the 
particle-hole representation. We studied systems of size 
L ranging from 2^" to 2^^. The system was allowed 
to evolve through 10^ — 10^ configurations before mea- 
surements were made. In order to speed up the algo- 
rithm to locate the site with the minimum /, we used a 
logarithmic-bin search procedure. Steady state averages 
were computed using ~ 10* configurations. While this 
number of configurations was averaged over to get an ac- 
curate estimate of the decay exponent, even an average 
over about 100 configurations indicates the presence of a 
strong density inhomogeneity clearly. 

In the untilted case, the height profile is an even func- 
tion (Fig. la) and '^{r) is an odd function (Fig. 10a) 
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FIG. 10. Density profiles in the (a) untilted (p = 0.5) 
and (b) tilted (p = 0.75) cases. Note that in (b), the value 
at the origin is larger than 0.75, indicating that the active 
site is more than nominally likely to have a particle, because 
of the bootstrap effect discussed in the text. While (a) is 
an odd function of r, (b) has no parity and the even and odd 
parts can be studied separately. The inset shows the odd part 
of the profile in the untilted (filled circles) and tilted (open 
circles) cases. Both decay as power laws though with differ- 
ing exponents. We used L = 16384 and averaged over 10® 
configurations. 

decaying asymptotically as a power law \r\~^ with 9 = 
0.90 ± 0.03 (Fig. 10, inset). In the tilted case as there 
is no r — > —r symmetry, \E'(r) does not have a definite 
parity (Fig. 10b). It is useful to separately analyse the 
even and odd functions 4'±(r) = (^'(r) ± *(-r))/2. The 
odd part decays as \l>_(r) ~ \r\~^- with 9- = 1.04±0.05 
(Fig. 10, inset). The even part *+(r) « -b{L)+a \r\~^+ 
where = 0.46 ± 0.05 and b{L) — » as the lattice size 
L ^ oo (Fig. 11). 

The density profiles of Fig. 10 correspond to the height 
patterns shown in Figs. 1 and 2. Qualitatively, the rea- 
son for this time-averaged structure of the height profile 
can be seen as follows. In the untilted case, on average, 
the active site is located at the peak (Fig. la) where fk 's 
which have not been sampled earlier are most likely to 
occur; such a region is thus more likely to contain small 
values of /. In the tilted case there is, in addition, a 
bootstrap effect at work. If p > 0.5, the active site is 
more likely to contain a particle. Given the dynamics, 
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FIG. 11. A log-log plot of the Fourier transform of the 
even part of the density profile in the tilted (p = 0.75) case, 
for the data appearing in Fig. 10. At low q the function de- 
cays as a power ~ g^*^ with ~ 0.54 implying that at large 
r , the behaviour is a power-law with decay r^'-^^'*-'. 

regions to the left are more often refreshed, making the 
active site likely to move in this direction and find itself 
amidst a particle cluster. The directionality of the active 
site motion is evident in Fig 6. This leads to densities 
higher than p on both sides of the active site (Fig. 10b) 
for p sufficiently different from 0.5. That particles are 
picked more often than would be expected on the basis 
of the nominal density is evident from the fact that the 
value of the density at the active-spot is larger than p. 
For the density pattern in Fig. 10b, which corresponds 
to p = 0.75, the value of the density at the origin is 
*(0) = 0.84. 

In the next section we develop an approximate theo- 
retical description for the exponents describing the decay 
of the pattern in terms of the exponent tt which describes 
the long-ranged hops of the active site. 



V. THE INTEGRAL EQUATION FOR PATTERN 
FORMATION 

As discussed in the previous section, a pattern in height 
gradients or densities is formed in the interface depinning 
model. It is related to the motion of the active site which 
rearranges the particles and holes in its wake. For exam- 
ple, if the active site were to remain stationary at a given 
site for a certain length of time, the dynamics is such that 
it would build up a pile-up of particles to the left of the 
site and a pile-up of holes to the right i. e. it would create 
a density shock around itself. In this section, we write 
down an integral equation which provides a description 
of the pattern in terms of p{l), the probability that the 
active site hops a length I in consecutive instants, as well 
as the local dynamics of interface readjustment. 
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If the pattern is centered at R{t) at time t, the dynam- 
ics causes two changes at the next instant: (i) A short- 
ranged readjustment of the interface changes the profile 
near R{t). The average density change at site R{t) + r 
is modelled through a density-increment function <i>(r) 
which is short-ranged, (ii) The active site jumps a dis- 
tance I = R{t+1) — R{t). Since the pattern is centered at 
the active site, the result of (i) followed by (ii) is that the 
average profile reproduces itself, except that it is centered 
at the shifted site R{t) + 1. Both effects are incorporated 
into the integral equation 

/+00 
[■i'{r-l) + <i>{r-l)]p{l)dl. (3) 
-oo 

This equation can be solved using Fourier transforms. 
In particular, the long-distance behaviour of the pattern 
^(r) is related to the decay of p{l) for large I, resulting in 
a scaling relation for the exponent 9 in terms of the ex- 
ponent TT. The details of the analysis are given in section 
V B below. 

However, this equation provides only an approximate 
description. To understand the nature of the approxi- 
mation made, we define a "Levy-flight model" , which is 
similar in spirit to the models of . As in these mod- 

els, there is no explicit quenched disorder, but the effect 
of disorder is modelled by a long-ranged jump probability 
distribution. For this model, we show that the integral 
equation (3) holds exactly. Further, we can gain consider- 
able insight into the mechanism of activity-centered pat- 
tern formation by studying the effect of different decays 
of the jump function p(l) within the Levy-flight model. 

A. Derivation of the Integral Equation for the 
Levy-fiight model 

The Levy-flight model is defined as follows. Consider 
a one-dimensional lattice with particles and holes. We 
assume that a jump probability distribution p{l) is spec- 
ified a priori: if at t = the active site is located at a 
site r, at the next instant, it can lie a distance I away 
{i.e. at I + r) with a probability p{l). Evidently there 
are no temporal correlations in this model, since at every 
instant, the jump length is chosen afresh from the distri- 
bution Spatial correlations in the active-site motion 
are however built in by hand since p{l) is given. Once a 
particle or hole is picked for update, the local rules are 
assumed to be the same as in the Extremal Bond Model 
i.e. the particle (or hole) exchanges position with the 
nearest hole (or particle) to the left (or right). 

We now derive an integral equation starting from the 
master equation for the Levy flight model. We first define 
a configuration i to be a set of integers {{ni{r)}, Ri) r — 
1...N where N is the size of the lattice and R is the cur- 
rent position of the active site. The variables n(r) can 



take the values or 1 depending on whether the site r is 
occupied by a particle or is empty and R can take any 
value between 1 and N. The total number of states is, 
therefore, Nt ~ N x ^ Cn^ where Np is the number 
of particles. In the usual manner, we characterise the 
steady state by a column vector \P). The entries of this 
column vector arc the probabilities Pi for the configura- 
tion i where i — 1....,Nt. To obtain the steady state, 
we need to solve the master equation d\P)/dt ~ W\P) 
where W is an Nt x Nt matrix connecting the different 
states. The dynamics that connects different states is 
the following: the particle (hole) at R is exchanged with 
the nearest hole (particle) to the left (right) as specified 
earlier. Subsequently the active site now hops from R 
to R + I with a probability p{l). The diagonal elements 
of the matrix are Wu = —1. The off-diagonal elements 
are given by Wij = p{l) if configuration i is connected 
to configuration j by an elementary update and a jump 
of the active site of length I and by Wij = otherwise . 
Since the probability p(Z) is normahsed, (X]p(0 = 1) f^e 
sum J2i foi' sach column of the matrix adds up to 0. 
This is a requirement for any stochastic matrix. 

Every configuration i can go to N other configurations. 
Each of these differ from i in the positions of the parti- 
cle and hole exchanged in the elementary update move, 
and also in the position of the active site. Similarly there 
are N configurations which feed into any configuration i. 
The construction of these configurations is very similar 
to that carried out in the case of the low-noise Toom- 
interface and related models |Q. However, unlike in 
that case, here the transition probabilities are not all the 
same. As a result the steady state here is very different 
from the product-measure steady state found in p9| ] and 
is difficult to characterise for general p{l). 

However we are able to characterise one aspect of the 
steady state by defining a quantity ^' in the following 
manner : 

Nt 

•^{r)^Y.P,n,{r + R,)-p (4) 

i=l 

That is, we add up the densities at a site r for each 
configuration that occurs in the steady state, having first 
shifted the origin separately in each configuration so that 
the site r for every state is one which lies r away from the 
active site. From this is subtracted the average density 
p. The above equation is equivalent to Eq. ^, with the 
time average in the latter being replaced by a weighted 
sum over configurations here. 

From examination of the W matrix, it is clear that in 
the case that p{l) is a constant independent of / (as in 
ordinary stochastic processes), all the P^'s are equal. As 
a result, from the definition, the ACP vanishes. How- 
ever, while this is a sufficient condition it is not neces- 
sary. There are instances (see Section VIII) when special 
symmetry considerations rule out any pattern formation. 
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Here 4>{r) is a short-ranged function (this is related to the 
'density-increment' function appearing in Eq. ^ which 
depends on the local update rules. It is defined by Eq. 
^ and is the difference between the two configurations i 
and j which are related by an update. This function is 
the same for this model as for the Bond model (Fig. 12) 
since it depends only on the local update rules. In the 
EBM, it is a short-ranged function whose range is deter- 
mined by the average length of particle and hole hops. 
Later in this section we will comment on circumstances 
in which this function can become long-ranged. 

Substituting Eq. ^ in Eq. H, we find that the right 
hand side can be rewritten as 

*M = E^^(0 T.^ri,{r + Mr " " P]P^ (7) 



Using the definition of ^'(r) again, wc finally obtain the 
integral equation |^ for 'l'(r) in terms of p{l), 



*(r) = ^p(/)[*(r - /) + $(r - /)] 



(8) 
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FIG. 12. The density-increment function $(r) in the EBM 
for (a) untilted and (b) tilted cases. Due to particle-hole sym- 
metry, the function is a short-ranged purely odd function in 
(a) while it has no specific parity in (b). In both cases how- 
ever, J ^(r)dr = from particle-hole conservation. We used 
L = 16384 and averaged over 16.10* configurations. 

However, the presence of an ACP definitely implies the 
presence of a nontrivial p{l); a pattern indicates correla- 
tions in the active-site motion. 

We will now derive Eq. ^ for the Levy-flight model. 
Note that the probability Pi = 'l2iPi^)^'j^i where 
the prime on the summation implies that only those con- 
figurations j are considered in which the active site is I 
sites away from its position in i. Further j should trans- 
form to i when the active-site is updated. This equation 
follows from the master equation. The sum is over all Vs. 

We can now substitute this in the RHS of Eq. ^ 



P 



(5) 



The density at site r in the configuration i is multiplied 
by the sum of the probabilities of those configurations 
that lead to it after a local update; a local exchange of a 
particle and hole, and a subsequent jump of the active- 
site of length /. Hence we can write 

n^{r + R,) ^ nj{r + R, - I) + (f)j{r + R, - I) (6) 



Here we have defined an averaged function (E>(r) — 
'^iPii'ii'"')- From particle conservation, it follows that 
J <^{r)dr = 0. For the symmetric case of half-filling, N/2 
particles and N /2 holes, the function $ is strictly an odd 
function. However this does not hold when the number 
of particles is not equal to the number of holes (Fig. 12 
b). 

Thus we have been able to show that the integral equa- 
tion (2) is valid for the Levy-flight model. From the na- 
ture of the model, it is clear that while spatial correla- 
tions in the active site motion are built in by hand, there 
are no temporal correlations in the length of subsequent 
hops of the active site. Therefore the integral equation is 
exact only in such a case. In the next section, however, 
we carry over some of the predictions of the integral equa- 
tion to the EBM and find that in some cases, it tallies 
quite well with numerical results. In Section VI, we dis- 
cuss briefiy how to generalise Eq. ^ to include temporal 
correlations such as are present in the EBM. 



B. Analysis of the integral equation 

We now investigate the predictions of Eq. ^ for ^'(r) 
in terms of a given active site hopping probability p{l). 
and the short-ranged readjustment function (j){r). We 
solve the equation using fourier transforms. Defining 
#((j) = e2'^'«'^*(r)(ir etc we find 



^(g)p(g) 



(9) 
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FIG. 13. The activity-centered pattern ^(r) in the Levy 
flight model for which p{l) decays rapidly as As evident 

in the figure, there is a near complete segregation of particles 
and holes about the active site. A system of size 1024 was 
considered and ~ 10^ configurations were averaged over. 



Since we are mainly interested in the q — > behaviour of 
Eq. H, we do not need the full functional form of ^{q) 
but only the leading order behaviour. 

Given a function p{l) the integral equation predicts 
a corresponding ^(r). The large-distance behaviour of 
the ACP thus depends on whether p{l) is short or long- 
ranged. We consider now three different cases for the 
function p{l) and solve for ^{q) using equation ^. We 
substantiate the predictions of the equation by numeri- 
cally simulating the levy flight model. 

Case 1: Consider first the case of an infinite-ranged 
p{l). The simplest case is when p{l) = jj where TV is the 
number of sites in the lattice. This case corresponds to 
usual stochastic processes. In this case p{q = 0) = 1 and 
p[q ^ 0) = 0. This implies that = for g ^ and 
is therefore a very short-ranged function in space. 

Case 2: Consider now the case when p{l) is a short- 
ranged function. We will consider the case when it is 
symmetric and hence for low q p{q) ^ q^ . Substitut- 
ing this in Eq. || we find that ^'(g) ~ sgn {q)/\q\. This 
implies that \E'(r) ~ sgn (r). The particles and holes sep- 
arate out completely and the active-site is located at the 
boundary between the two. This is easy to understand if 
we consider the limiting case when the active site is to- 
tally stationary. In this case, the action of the dynamics 
is to move all the particles from the right of the active 
site to its left. Eventually, this leads to a total separa- 
tion of particles and holes. This picture is modified only 
slightly when the active site executes a localised motion 
about any lattice site, and hence for a short-ranged p{l) 
(Fig. 13). 

Case 3: We now come to the case of interest for the 
interface depinning model i.e. when p{l) decays as a 
slow power law. If p+{l) = , Eq. ^ predicts that 

\l/(r) ^ sgn{r)\rf- . The exponents tt and 9- are related 



by the scaling relation 9- + = 3. This is to be com- 
pared with the numerical estimate 0_ -t-Tr^ = 3.15 for the 
Extremal Bond Model in the untilted case. Consider now 
the case when none of the functions p{l), vE'(r) and <I>(r) 
has a definite parity. This is relevant for the tilted inter- 
face in the Sneppcn model. Since $(r) is short-ranged, 
$(q) ~ i(j)iq + (j)2q^ as q ^ 0. There is no (po term, as the 
elementary step of hopping a particle or hole conserves 
particle number, implying J ^(r)dr ~ 0. The q ^ be- 
haviour of p{q) is determined by the asymptotic power 
law decays of the even and odd parts p±{l) as \l\ — > oo. 
Thus we have p~+{q) « 1 — We might have 

expected p'-iq) ~ Bq + C sgn{q) \q\'^-~^, but in fact 
the mean velocity / lP{l)dl of the active site vanishes (as 
mentioned in Section III) implying B = 0. Thus the inte- 
gral equation predicts that to leading order, both ^'+(r) 
and 5'_(r) decay as powers ~ |r|^^±, with 



and 



2tt+ — 7r_ = 3 



+TT+ 



The prediction ir^ + 6^ =3 compares quite well with 
the numerically determined values 3.04 for ^'_(r) in the 
tilted case. For ^+(r), however, the numerically deter- 
mined value of 9+ + 27r+ — 7r_(~ 1.97) deviates substan- 
tially from the predicted value 3. The likely reason be- 
hind the discrepancy is explained in the next section. 



C. Drawbacks of the Approximation 

The integral equation (2) is exact for the Levy flight 
model, but only approximate for the EBM of interface 
depinning. Here we briefly run over the nature of the ap- 
proximations made, and possible directions for improve- 
ment. 

One sort of approximation is the neglect of correlations 
between lengths of successive jumps; as we saw in Section 
III, these correlations are very marked. The extension of 
the integral equation to include such correlations is dis- 
cussed in Section VI B. 

As noted in Case 3 of Section VB for tilted interfaces 
{p / 0.5), the pattern is not well represented by Eq. (2) 
even qualitatively. The reason for this is that the p ^ 0.5 
pattern in the EBM is formed due to a feedback effect 
which is missing in the Levy flight model. If p > 0.5, then 
within EBM dynamics, particles (or holes if p < 0.5) are 
picked more often than just p.Ntry times in Ntry tries. 
This condition is not incorporated in Eq. ^ at all. To 
test how important this effect is, we simulated the Levy 
flight model with the further constraint that if a site cho- 
sen for growth does not contain a particle, it is discarded 
and the search is continued till a site containing a par- 
ticle is found. This leads to a pattern that much more 
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FIG. 14. The activity-centered pattern 'l'(r) for a Levy 
Flight Model with the further constraint that only sites con- 
taining particles are chosen. The pattern generated with these 
rules resembles 'i'+{r) in the EBM more closely than if the 
rules were implemented without the constraint. This fits in 
with our conjecture as to why the Integral equation does not 
describe the even part of the pattern. 

resembles that in the p ^ 0.5 EBM, in that the even part 
of the pattern is much more prominent than the odd part 
(see Fig. 14). 

Another point to note is that we have always taken 
the density-increment function $ to be short-ranged. In- 
sofar as we are interested in the EBM, this is certainly 
so. However, in the Levy-flight model, any power larger 
than 3, results in a pattern like the one shown in Fig. 13, 
with a segregation of particles and holes and not one de- 
caying as a power as predicted by the integral equation. 
The reason for this is that in this case, the mean squared 
distance (i?^(t)) covered by the active site in time t is fi- 
nite. Hence for a large system, this is like a short-ranged 
motion and leads to particle hole segregation. In the pro- 
cess, as particle clusters build up, the density-increment 
function is no longer short-ranged. 

Lastly, while in the EBM, the probability distribution 
of jumps p{l) arises naturally from the dynamics, in the 
Levy-flight model we put in the spatial correlation in 
the active site motion by hand. In the following section, 
we try and remedy this point by modifying the integral 
equation by writing a set of coupled equations. 



largest probability of being visited again since it is most 
likely to have the least random number. The more often 
a region is visited, the more the random numbers associ- 
ated with sites in that region are refreshed and hence the 
more likely it is that the minimum random number lies in 
that region. Since the dynamics creates a particle excess 
on the left of the active site and a corresponding paucity 
on its right, a density shock is built up in that region. 
The extent of the shock is directly proportional to the 
length of time spent by the active spot in that region. 
The larger the shock, the larger the density difference, 
and hence the larger the derivative of 5*. Hence the ef- 
fect of a region being visited a large number of times is 
precisely that /dr is large. But as explained above, 
the larger the number of times a site is visited, the larger 
is the probability of its being visited again. Thus we con- 
jecture that the probability of visiting a region a distance 
I away from the active site is proportional to the mag- 
nitude of d"^ / dr in that region, and write the following 
equations: 



*(r) = / dlp{l)[^{r + $(r - I)] 



and 



dr 



(10) 



(11) 



where A/" = / |d5'/dr| is the appropriate normalising fac- 
tor. 

Equation |l^ is the integral equation that we derived 
in the last section. Equation 1^ now defines p{l) in 
terms of the pattern "if. The details of the dynamics of 
particle-hole exchange in this model make this a reason- 
able premise as explained above. 

Let us compare the results of the coupled equations 
with those for the untilted interface. If ^(r) ^ \r\~^ and 
p{l) ~ 1^1"'^, we saw in section VB that Eq. |l^ im- 
plies relation tt -I- 6* = 3. Further, Eq. |ll| predicts the 
relation 6 + 1 = tt. The solution is tt = 2 and 9=1. 
This compares quite well with the value of the exponents 
TT = 2.25 ± 0.05 and 6* = 0.9 ± 0.03 found numerically for 
the EBM. 



VI. TEMPORAL CORRELATIONS IN THE 
INTERFACE DEPINNING MODEL 



D. The Coupled Equations 



A. Time dependence of the pattern 



A drawback of the Levy flight model is that the func- 
tion p{l) has to be supplied from the outside. We try 
to rectify this by writing a set of coupled equations that 
predict ^(r) and p{l) in terms of each other. 

In order to do this, we need to understand the proba- 
bility that a certain portion of the pattern is visited by 
the active site. In the EBM, the site visited last has the 



The pattern, as defined earlier, is obtained by a 
straight time average over all configurations, taking care 
to shift the origin appropriately at every instant. It is 
interesting to ask how the pattern changes if we con- 
sider only those configurations that occur after very large 
jumps of the active site (or equivalently, only those con- 
figurations for which the minimum random number is 



13 



0.9 



0.8 

VF(r) 

0.7 

















"""^ |-^iTtDiH twmti 







-200 200 

r 



FIG. 15. Evidence for the dependence of the AGP on the 
length of the previous jump of the active site, ^'(r) is shown 
with time-averaging done over (1) all configurations (open 
circles) (2) only configurations resulting after a jump in the 
active site of magnitude 50 (filled triangles) (3) only configu- 
rations resulting after a jump of length 100 (open triangles) 
and (4) only configurations resulting after a jump of length 
300 (filled circles). As is evident, the pattern gets more and 
more squashed till it barely exists in (4). A system of size 
1024 was used and ~ 10^ configurations were averaged over. 

very close to the largest possible). A strong dependence 
of the pattern on the configurations averaged over is ev- 
ident in Fig. 15. The pattern gets more squashed as 
the jumps in the active site leading to the configuration 
get larger. This implies that the pattern ceases to exist 
at stoppers and builds up again as the interface pierces 
through. The actual dynamics of pattern collapse and 
build-up is an interesting subject for further study. 

It should be recalled that long jumps occur much more 
infrequently than short jumps (Fig. 7); the probability of 
occurrence of a stopper is thus very low. The time aver- 
age in the definition of the pattern (Eq. |^) is dominated 
by configurations between stoppers (when the active site 
is moving within the loops of the directed percolation 
network) , rather than those at stoppers (when the active 
site is on the backbone of the network). 

B. A Hierarchy of Integral Equations 

The integral equation 3 can be modified to include cor- 
relations in time. To do this, we enlarge the definition 
of a configuration i to include the active site at the pre- 
vious instant as well. Now the off-diagonal elements of 
the transition matrix elements Wij — p{l'\l) if configu- 
ration j has resulted as a consequence of a jump of the 
active site of length I and is connected to configuration 
i by an elementary update and a jump of the active site 
of length I'. The function p{l'\l) is just the conditional 
jump probability already introduced in Section III. 

Following the same procedure as before, we get an 



equation for Vf//, the pattern resulting from an active 
site hop of length I' : 

*i'(^) = - n + Mr - 0] • (12) 

Keeping correlations upto one time step back gives the 
pattern a non-trivial dependence on the jump length 
The integral equation || gives only a trivial dependence 
of ^{r) on I'. 

This procedure can be further generalised by going 
back one more step in time and keeping the location of 
the active site two instants back. 

This leads to an equation of the following sort: 

W = Y.P(^\^M'^h,hir -l) + $;„i,(r - I)] 
h 

where p{l\li\l2) is the conditional probability that the ac- 
tive site hops a length I given that at the previous instant 
it hopped a distance h and in the instant before that, a 
length I2. Similarly, '^ii,i2{'r') is the pattern formed when 
averaged over configurations that result after two consec- 
utive jumps of h and I2 respectively. 
Here J^h "^Uiir) */('')■ 

Keeping the time sequence of jumps leads to an infinite 
hierarchy of equations. 

= E E?'(^'io[*''('' - + *''(^ - 0] (13) 

^lir) = E?'(^l^il^2)[«';„i.(r ~ /) + <i>,„i,(r ~ 0] (14) 

h,l2 

= (15) 

The integral equation ^ corresponds to curtailing this 
hierarchy at the first step by assuming that 5';(r — I') = 
^{r - I) and p{l'\l) = p{l' . 



VII. TWO-POINT CORRELATION FUNCTIONS 
IN THE INTERFACE DEPINNING MODEL 

The activity-centred pattern ^(r) defined in Eq. 1 is 
a one-point correlation function, defined with respect to 
an erratically-moving origin. As a result, it has a strong 
effect on the customary, space-time averaged two-point 
correlation function 

C(Ar) = {(n(r')n(r' + Ar))} - 

Here r' is a fixed site on the lattice, (...) stands for a time 
average and {...} stands for an average over all sites r' . 

Numerical results for C(Ar) for a tilted interface show 
that it saturates at a value C""* which decreases with 
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Ar 



FIG. 16. Two-point density- density correlation function 
for L — 4096 (open circles) and for L — 16384 (squares). The 
saturation value is reduced strongly (filled circles) on sub- 
tracting the contribution of the AGP. This lends support to 
the idea that the one-point correlation function defined by 
the AGP enters the definition of the connected part of the 
two-point correlation function. The data was obtained by 
averaging over 10® configurations. 



as Ar — > oo. Now consider averaging the function 
r(r, Ar) over r. On performing a space average over the 
right hand side of Eq. |l6|, we note that 

{(Sn{r+R{t))Sn{r+Ar+R{t))}} = {{6n{r)Sn{r + Ar))}. 

Thus this implies that {(n(r)n(r -I- Ar))} — 
{(/9 -t- ^'(r))(p -f vE'(r -I- Ar))} approaches zero as Ar 
oo. This predicts the saturation value of the correlation 
function 

C'"' = {{p + *(r))(p + *(r + Ar))} - p\ 

To test this, we subtracted this estimate of the saturation 
value from C(Ar) and found that the saturation effect is 
in fact suppressed strongly (Fig. 16), supporting our in- 
terpretation. It is possible that another slightly different 
definition of the pattern would eliminate the slight shoul- 
der which remains in Fig. 16, after subtraction of C*"*. 

Usually, if the two-point correlation function saturates 
as the separation between the two points is increased, the 
saturation value is associated with a nonzero value of the 
space-fixed average {n{r')) — p. The unusual aspect here 
is that there is saturation even though {n(r')) = p. 

We studied the manner in which C(Ar) — C^"* ap- 
proaches zero as Ar oo. In the tilted case, we found 
that the correlation function decays exponentially, by 
studying a model where the rules are the same except 
that only particles are picked. This corresponds to a case 
of extreme tilt and is similar to the case studied by . 
In the untilted case, C(Ar) — C""* decays as a power law 
~ r^'^ with K ~ 0.6. 



increasing size L ^ (Fig 16). This happens because sub- 
tracting out the quantity p^ from the correlation func- 
tion is not correct, since the presence of the ACP causes 
a density inhomogeneity in the medium which cannot be 
accounted for by subtracting out a constant quantity. To 
account for the presence of the pattern we consider the 
correlation function 

r(r, Ar) = (5n{r + R{t))Sn{r + Ar + R{t))) (16) 

where R{t) is the location of the active site , r is the dis- 
tance from the active site and dn{r+R{t)) = n{r+R{t)) — 
p — vlf(r) is the fluctuation around the average ACP. A 
reasonable expectation is that the (5n's are independent 
for large separations Ar, i.e. 

r(r, Ar) = {6n{r + R{t)) {Sn{r + R{t) + Ar)) ~> 



VIII. PATTERNS IN OTHER MODELS 

In the Extremal Bond Model of interface depinning, 
we have described the activity-centered pattern in height 
gradients which forms as a result of correlated motion of 
the active site. However, this is not the only sort of pat- 
tern that is formed. There is pattern formation also in 
the value of the average random number / at a site, as a 
function of the distance from the active site. In analogy 
with Eq. 1^, we define this pattern in random numbers as 

^fir) = {fir + R{t)))-{{f)}. 

Here /(r -I- R{t)) is the random number at a distance r 
from the active-site and the time average (...) and space 



'^At much larger values of r, proportional to the system size 
L, C{r) changes sign, as a result of the particle conservation 
sum rule J C{r)dr = 0. 
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FIG. 17. The pattern in random numbers for the EBM 
(the untilted case) in (a) real space (b) Fourier space. The 
Fourier transformed function appears to diverge as a power 
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-(1-0) 



at 



large r for the function in real space. The above data is for a 
system of size L = 8192, averaging over 10^ configurations. 

average {...} are performed over configurations in the 
steady state, as before. 

Moreover, patterns are found in other extremal mod- 
els as well. Figures 17, 18 and 19 show the /-patterns 
in the EBM, the Bak-Sneppen model of biological evo- 
lution [p8[ and the Zaitsev model of low-temperature 
creep |37]] respectively. Numerically, it is difficult to di- 
rectly extract the manner in which the patterns shown 
in Figs. 17-19 approach their asymptotic values, as fits 
to power law decays are very sensitive to the assumed 
saturation value. We avoided this problem by studying 
the Fourier transforms of the functions, as the saturation 
value influences only the single Fourier mode at q = 0. 
In all three cases, the Fourier transforms show evidence 
of power law behaviour as g ^ 0, implying power law ap- 
proaches of the /-patterns to their respective saturation 
values in real space. However, we have not developed an 
analytical description of the /-patterns in any of these 
models. 

Pattern formation thus seems to be generic to extremal 
models. However there are instances when considerations 
of symmetry rule out the formation of a pattern. This 
is exemplified in the following model, similar in spirit to 
that considered in |5^. The rules of the dynamics are 
the same as in the EBM, random numbers are assigned to 



FIG. 18. The pattern in random numbers for the 
Bak-Sneppen model for evolution in (a) real space (b) Fourier 
space. The Fourier transformed function appears to diverge 
as a power g"** with (f> ~ 0.24 implying a decay exponent 
(1 — (j)) for the function in real space. This pattern could sig- 
nify, for instance, the average fitness of a species as a function 
of the distance from the currently mutating one. The above 
data is for a system of size L = 8192 and 10^ configurations 
have been averaged over. 

every site and the minimum is picked, except that there 
is no net current of particles as there was in the EBM. 
If the site picked is occupied by a particle, the particle 
exchanges place with the first hole to its left, and if the 
site is occupied by a hole, the hole too exchanges place 
with the first particle to its left. At half filling, there is no 
net current. This model is the most symmetric of those 
considered so far and there is no density pattern formed. 
However, as in the Bak-Sneppen and the Zaitsev models, 
there is a nontrivial pattern in random numbers in this 
model. 

We emphasize that the feature of the dynamics which is 
responsible for activity-centered pattern formation is the 
existence of correlations in the motion of the active site. 
Extremal models constitute just one class in which there 
are such correlations. An example of another such class 
is models of certain types of reaction-diffusion systems, 
where the activity is quite constrained and correlated. 
In another physical context, it would seem that coherent 
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FIG. 19. The pattern in random numbers for the Zait- 
sev model in (a) real space (b) Fourier space. The Fourier 
transformed function appears to diverge as a power g"*^ with 
(j) ~ 0.26 implying a decay (1 — phi) for the function in real 
space. The data shown is for a system of size L — 2048, with 
an average over 10^ configurations. 

structures which form in turbulent flows [ pT| may well be 
described by activity-centered patterns. In our definition 
of the pattern, we take an average over all times, keeping 
track of the moving structure. This is to be contrasted 
with previously used methods to identify such moving 
structures, based on the notion of conditional sampling, 
namely, averaging only over those time zones in which the 
activity is at a particular space-time location |5l| ]. Our 
definition takes configurations at all times into account, 
but requires a shift of the origin at every instant. 



IX. SUMMARY AND CONCLUDING REMARKS 

We have introduced a variant of the Sneppen model 
of interface depinning — the Extremal Bond Model — 
and have studied the effect of the dynamics of the growth 
process on the shape of the interface. Our principal re- 
sult is the observation that there is a nontrivial structure 
which forms in the interface, and which moves along with 
the active site. A simple time average of the height gra- 
dients, measured in a frame of reference which moves 



with the active site, defines the activity-centered pattern 
which serves to quantify the structure. Our numerical 
study shows that the pattern has a tail which decays as 
a power law at large distances. 

An understanding of the mechanism underlying 
activity-centered pattern formation was obtained by 
writing an integral equation which relates the pattern to 
the probability distribution of active-site jump lengths. 
The integral equation could be derived by writing an ex- 
tended master equation for the Levy flight model, making 
it clear that the equation is exact only when there are no 
temporal correlations between successive jump lengths. 
In the extremal bond model, however, such correlations 
are strong. We have shown that a correct description of 
the pattern, involving temporal correlations, necessitates 
keeping an infinite hierarchy of equations. Terminating 
this hierarchy at the very first step results in our inte- 
gral equation. It would be of interest to understand how 
much one can better this description by keeping more 
steps in the hierarchy. 

The activity-centered pattern is a one-point function, 
and so enters the definition of two-point correlation func- 
tions. The physical point is that the density inhomogene- 
ity caused by the pattern must be taken into account, 
by subtracting the relevant quantity from the density- 
density correlation function. If this is not done, and the 
square of the the nominal density p is subtracted instead, 
the correlation function can exhibit a rather unusual sort 
of finite size effect. This is an interesting point since the 
space-fixed time-averaged density at a site is p, and it 
is only the density as defined in Eq. ^ with respect to 
the moving active site which is different from the nom- 
inal density. Yet the space-fixed two-point function is 
affected by this "hidden" pattern. 

The presence of the pattern clearly points to a non- 
homogeneity in the interface: the region around the ac- 
tive site looks very different on average from the region 
far away from it. For instance, we expect there to be a 
larger length of interface in a region of fixed size x around 
the active site, than in a region opposite it. We mon- 
itored mean squared fluctuations of the height around 
the instantaneous average, in regions around and oppo- 
site the active site, and found a pronounced difference 
(factor ~ 2, for both tilted and untilted interfaces with 
X = 256, L = 4096). This effect is smaller at stoppers, 
in keeping with our finding that the pattern itself is sup- 
pressed there. This excess length of interface associated 
with the activity-centered pattern may provide a useful 
way to identify the active region in experiment. 

Finally, it was pointed out that activity-centered pat- 
tern formation may occur in a wide variety of other phys- 
ical contexts, ranging from low-temperature creep of dis- 
locations to structures in turbulent flows. We have pre- 
sented numerical evidence for this sort of pattern forma- 
tion in a number of other extremal models. But more 
generally we expect activity-centered patterns to form 
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whenever there are strong correlations between succes- 
sive locations of the active site. 
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